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^ . Abstract 

' In order to investigate to what extent is the low-lying behavior of even-even nuclei dependent on 

, particular nucleon-nucleon interactions, we consider systems of bosons where these interactions are 

taken as gaussian random numbers with equal likelyhood of being attractive or repulsive. We find a 
' statistical dominance of L — ground states and other correlations, which we analyze in terms of a 

, mean field approach. 

Para investigar hasta que punto el comportamiento de los micleos esta determinado por interac- 
ciones particulares entre los nucleones, consideramos sistemas de bosones con interacciones gaussianas 
aleatorias, donde dichas interacciones tienen la misma probabilidad de ser atractivas o repulsivas. En- 
contramos un dominio estadfstico de estados base L — y otras correlaciones, que analizamos en 
terminos de un metodo de campo medio. 



PACS numbers: 05.30. Jp, 21.60.Ev, 21.60.Fw, 24.60. Lz 



1 Introduction 



Whereas conventional wisdom ascribes the observed properties of low-lying nuclear systems to a particu- 
lar form of the nucleon-nucleon interaction, recent investigations have found that some of these features 
may represent a general and robust property of many-body systems 0-0. Random matrices have been 
used in the past in order to study generic spectral properties, but the emphasis was usually centered on 
global characteristics, such as level-spacing distributions in highly excited nuclei or the average proper- 
ties of small metallic particles or quantum dots ||. This approach had its origin in the observation by 
E.P. Wigner that . . the Hamiltonian which governs the behavior of a complicated system is a random 
symmetric matrix, with no particular properties except for its symmetric nature '. At high energy, nuclear 
dynamics is thus assumed to have lost track of any correlations. In the case of other systems, random fluc- 
tuations in shapes and/or impurities can be reasonably modeled in this form and their average properties 
reproduced. 

A general conclusion of these studies is that the specific mechanism (interactions) seem to be irrelevant 
in the determination of global features. 

However, the random-matrix analyses of nuclear distributions typically involve states with the same 
quantum numbers, such as angular momentum, parity or isospin, and little attention was given to cor- 
relations among different symmetries. Although the so-called two-body random ensemble (TBRE) was 
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defined long time ago and some of its properties investigated p, |i0| , again attention was focused mostly 
on its global features and not so much on low-lying patterns H]7 This kind of ensemble is associated with 
a random selection of the two-body matrix elements, while the n-body Hamiltonian matrices are gener- 
ated with the usual coupling procedures. On the other hand, it is known that nuclei display very regular 
spectral properties. For example, in Refs. p"l| , |l2] ] an analysis of energy systematics of medium and heavy 
even-even nuclei suggests a classification in terms of seniority, anharmonic vibrator and rotor regions. 
Plots of the excitation energies of the yrast states with L p = 4 + against L p = 2 + show characteristic 
slopes for each region. This and other correlations have been shown to be robust features of low-energy 
nuclear behavior which signal the emergence of order and collectivity. In each case this behavior has 
been shown to arise from particular nucleon-nucleon interactions, such as an attractive pairing force in 
semimagic nuclei and the strongly attractive neutron-proton quadrupole-quadrupole interaction for de- 
formed nuclei. It came as a surprise, therefore, that recent studies for even-even nuclei, in a shell model 
with randomly distributed two-body interactions (TBRE), displayed a marked statistical preference for 
L p = + ground states, energy gaps, and other indicators of ordered behavior Q. These results sparked 
a large number of other numerical and theoretical investigations, in both fermion ||, ||, [7], [ll], [l4| and 
boson systems [| 0, |l5[ 0, in an attempt to explain this remarkable and unexpected results. 

In this contribution we summarize some of the most notable findings in the context of the Interact- 
ing Boson Model (IBM) and present a possible explanation in terms of a mean-field analysis. In the 
conclusions we indicate some of the many remaining open questions and future research. 



2 The interacting boson model 

The IBM describes low-lying collective excitations in nuclei in terms of a system of N interacting monopole 
and quadrupole bosons [ jl8| |. An analysis of the IBM with random one- and two-body interactions gave 
rise to further statistical evidence of spontaneous emergence of order, manifested by a large percentage 
of ground states with L p = + and clear vibrational and rotational correlations j|, |). Fig. [I] shows that 
there is a predominance (63.4 %) of L p = + ground states; in 13.8 % of the cases the ground state has 
L p = 2 + , and in 16.7 % it has the maximum value of the angular momentum L p — 32 + . For the cases 
with a L p = + ground state, we show in Fig. ^ the probability distribution of the energy ratio 

There are two very pronounced peaks, right at the vibrational value of R = 2 and at the rotational 
value of R — 10/3, a clear indication of the occurrence of vibrational and rotational structure. This is 
confirmed by a simultaneous study of the quadrupole transitions between the levels . 

These are surprising results in the sense that, according to the conventional ideas in the field, the 
occurrence of L = ground states and the existence of vibrational and rotational bands are due to very 
specific forms of the interactions. The study of the IBM Hamiltonian with random one- and two-body 
interactions seems to indicate that this may not be the entire story. However, the above results were 
obtained from numerical studies. It would be very interesting to gain a better understanding as to why 
this happens. What is the origin of the regular features which arise from random (both in sign and 
size) interactions? In this respect, there is a relevant quote by E.P. Wigner (as communicated to us by 
M. Moshinsky): 'I am happy to learn that the computer understands the problem, but I would like to 
understand it too'. 

In the next sections, we report on an attempt in this direction by considering the vibron model, 
which has the same qualitative features as the IBM, namely vibrational and rotational spectra, but has 
a much simpler mathematical structure, since there are no coefficients of fractional parentage (cfp's). All 
many-body matrix elements can be derived in closed analytic form [Q9l E0| . 
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3 The vibron model 



The vibron model is an interacting boson model designed to describe the relative motion in two-body 
problems, e.g. diatomic molecules [^J, nuclear clusters |^2j and mesons Its building blocks are a 

dipole boson p^ with L p = 1~ and a scalar boson with L p = + . The total number of bosons N is 
conserved by the vibron Hamiltonian. Here we only consider to one- and two-body interactions 
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where Hi contains the boson energies 



(2) 
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and H 2 all possible two-body interactions 
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and p m = (— l) 1_m p_ m . We have scaled Hi by TV and by N(N — 1) in order to remove the N 
dependence of the matrix elements. 

We first carry out a numerical study of the vibron model using random one- and two-body interactions. 
The seven parameters of the Hamiltonian altogether denoted by 



(•r) 



(e s , e p: u ,ui, Co, Co., v ) , 



are taken as independent random numbers on a Gaussian distribution 

P( Xl ) = e-^/^/aV^ , 



(5) 



(G) 



with zero mean and width a. In Fig. ^ we show the percentages of L = 0, L = 1 and L = N ground 
states as a function of the total number of vibrons N. Just as for the IBM, the vibron model shows a 
dominance of L — ground states. For even values of N the ground state has L = in ~ 71 % of the 
cases, and for odd values in ~ 55 % of the cases. Similarly, the percentage of ground states with L = 1 
shows an oscillation between ~ 1 % for even values of N and ~ 18 % for odd values. In ~ 24 % of the 
cases the ground state has the maximum value of the angular momentum L = N . 

For the cases with a L p = + ground state, we show in Figs. |J and |j] the probability distribution of 
the energy ratio 



R = 



E 2 + 



E ot 



E X7 



(7) 



In the vibrational limit the spectrum is that of a three-dimensional harmonic oscillator with R = 2, and 
in the rotational limit that of a three-dimensional Morse oscillator with R = 3. Both for even and odd 
values of N we see two clear peaks, one at the vibrational value of R = 2 and one at the rotational value 
of R = 3. Moreover, for even values of N there is a peak at R = 0, which is absent for odd values, as we 
shall prove below. 
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This numerical study confirms that the vibron model exhibits the same regular features as does 
the IBM, although there are some differences as well. Now the question is, how can we understand 
these properties in an analytic and more intuitive way? In a recent study the tridiagonal form of the 
Hamiltonian matrix of the vibron model was used to establish a connection with random polynomials 
]Tg| , |l7| ]. However, in general the Hamiltonian matrix is not of this form, and one has to look for alternative 
methods. 

In the next section we show that a mean-field study of the vibron model with random interactions 
can account for all features that were discussed above. 



4 Mean-field analysis 

The connection between the vibron model, potential energy surfaces, equilibrium configurations, shapes, 
etc., can be studied by means of coherent states. The coherent state for the vibron model can be written 
as a condensate of a deformed boson which is a superposition of a scalar and a dipole boson pit] 

\N,a) = -L= (y/l-a* S i+apl) N |0) , (8) 



with < a < 1. The potential energy surface is then given by the expectation value of the Hamiltonian 
in the coherent state 

En (a) = (N,a\ H \N,a) = a 4 a i + a 2 a 2 + a , (9) 
where the coefficients ai are linear combinations of the parameters of the Hamiltonian 

1 111 

a 4 = r-x = -u + ui + -c + -C2 + —f=v 0l 

1 

a2 = s ■ x = —e s + e p - it - u\ —vo , 

V6 

1 , N 

ao = e s + -jUo . (10) 

For random interactions, we expect the trial wave function of Eq. (||) and the energy surface of Eq. @ 
to provide information on the distribution of shapes that the model can acquire. The equilibrium con- 
figuration is characterized by the value of a = ao for which the energy surface En (a) has its minimum 
value. For a given Hamiltonian, the value of ao depends on the coefficients a\ and a-i- The distribution 
of shapes for an ensemble of Hamiltonians depends on the joint probability distribution of the coefficients 
0,4 and a 2 



P{ai,a,2) = J Hi dxi P(xi) 6(a4 — r ■ x) 5(ci2 — s ■ x) 
1 



2tt det M 

where P(xi) is given by Eq.dq) and 



cxp 



I(a 4 J) 



(11) 



r ■ s s ■ s 



M - ( I'l I'i ) • (12) 



The vectors r and s are defined in Eq. ( ^0| ) . In this approximation, we find that there are only three 
possible equilibrium configurations: 
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ao = 0: s-boson condensate. This corresponds to a spherical shape whose probability can be 
obtained by integrating P(a4, 02) over the appropriate range Si 



Pi = dci4 da-i P(ci4, 02) 

J Si 

1 f\s-s + f-s\ 

= — tt + 2 arctan ■ — 

4tt [ V VdetM 

This configuration has spherical symmetry and hence can only have L = 



= 0.396 . (13) 




• < a < 1: deformed condensate. The probability for the occurrence of a deformed shape can be 
calculated in a similar way and is given by 

/ o„ I am- S r \ 

= 0.216 . (14) 

A deformed shape corresponds to a rotational band with angular momenta L = 0, 1, . . . , N. 

• «o = 1- The probability for finding the third solution, a p-boson condensate, is given by 

P 3 = I-P1-P2 = 0.388 . 

The angular momentum content of a p-boson condensate is L — N, N — 2, . . . , 1 for N odd or 
L = N, N - 2, . . . , for N even. 

The probability that the ground state has L — can be estimated by assuming a complete decoupling of 
vibrations and rotations. For N even this gives ^(1 + Pi) or 69.8 %, and for N odd we have Pi + \Pi 
or 50.4 %, in good agreement with the numerical values. In Table |l| we show the probabilities and 
corresponding percentages for the ground states with angular momentum L — 0, 1 and N . 

A better estimate can be obtained by evaluating the moment of inertia. We adopt the Thouless-Valatin 
prescription, which leads to the formula 

2Na ° (15) 



4 (^-!) [575«o(l-a2)-i(co-c 2 )a2 
The moment of inertia thus depends in a complicated way on the parameters in the Hamiltonian, both 



explicitly as seen in the denominator of Eq. (15) and implicitly through ao, which determines the equilib- 
rium configuration. Formally, the probability to find a L — ground state is determined by the integral 
of the joint probability distribution P(fl4, 02, T) over the appropriate range. The results are summarized 
in Table ||. The spherical shape correponds to an energy ratio of R = 2, and the deformed shape to 
R = 3. This explains the vibrational and rotational peaks observed in Fig. [I] and [|. The occurrence of 
a peak at R = for even values of N is related to the p-boson condensate solution, which has angular 
momenta L = 0, 2, . . . , TV. The first excited L = 1 state belongs to a different band and has a higher 
excitation energy. For odd values of N the p-boson condensate has no state with L = 0, and hence the 
peak at R — is absent. 

In Fig. ^ we show the percentages of L = 0, L = 1 and L = N ground states, as calculated in the 
mean-field analysis. A comparison with the exact results of Fig. || shows that the mean-field results are in 
excellent agreement with the exact ones. The oscillations in the percentages of L = and L — 1 ground 
states are entirely due to the angular momentum decomposition of the p-boson condensate, whereas the 
percentage of ground states with the maximum value of the angular momentum L — N is a constant and 
does not depend on N. 
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5 Summary and conclusions 



In this contribution, we have studied the properties of low-lying collective levels in the interacting boson 
model and the vibron model with random interactions. In particular, we discussed the dominance of 
L = ground states and the occurrence of vibrational and rotational band structures. In both models 
these features represent general and robust properties, and are a consequence of the many-body dynamics. 

We have studied the origin of these regular features in more detail in the context of the vibron model. 
In a mean-field analysis of the vibron model with random interactions it was found that we can account 
for all observed features, both qualitatively and quantitavely. The occurrence of L = ground states 
can be related to the different equilibrium configurations of the vibron model. There are three different 
equilibrium configurations: a spherical shape (40 %), a deformed shape (20 %) and a condensate of 
dipole bosons (40 %). Since the spherical shape only has L = 0, and the deformed shape and the p- 
boson condensate with N even in about half the number of cases, one finds a L = ground state in 
approximately 70 % of the cases for TV even and 50 % for N odd. The spherical shape gives rise to the 
occurrence of vibrational structure, and the deformed shape to rotational bands. 

The advantage of the mean-field method is that it can be applied directly to the IBM as well. The use 
of coherent states circumvents the use of coefficients of fractional parentage, and allows one to associate 
different regions of the parameter space with geometric shapes. 

There are many open questions in this area of research. Although the results in the case of the collective 
boson models can be explained by the allowed shape distributions, there is still no clear explanation for 
the dominance of L = ground states in fermion systems in the context of the shell model. For example, 
a recent analysis for single j-shell random Hamiltonians displays strong oscillations for the L — ground 



state percentage as a function of j 14 1, reminiscent of the variations found for the IBM and the vibron 
model as a function of the boson number N || . The role of the approximate conservation of seniority 
should be further investigated. Another outstanding question is whether an appropriate truncation of 
the shell model space can lead to vibrational and rotational correlations of the kind observed in the IBM. 
Can shape be more precisely defined, for example, in a restricted fermion-paired Hilbert space? These 
and other studies are in progress. 
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Table 1: Probabilities and percentages for ground states with angular momentum L in the decoupling 
limit 



N even N odd 

Ground State Prob. Perc. Prob. Perc. 

L = 1(1 + Pi) 69.8% Pi + \Pi 50.4% 

L=l 0.0% i(l-Pi-P 2 ) 19.4% 

L = N 5(1 -Pi) 30.2% 5(1 -Pi) 30.2% 



Table 2: Percentages of L = ground states, calculated in the decoupling limit (mf) and with the moment 
of inertia (J). The exact results are based on 100,000 runs with N = 20 and N = 21. 
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N odd 




Shape 
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mf 
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39.6 % 


39.5 % 


39.6 % 


39.5 % 
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< a < 1 


10.8 % 


13.8 % 


10.8 % 


13.8 % 


3 


"0 = 1 


19.4 % 


17.9 % 


0.0 % 


0.0 % 





Total 


69.8 % 


71.2 % 


50.4 % 


53.3 % 




Exact 




71.2 % 




54.7 % 
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Figure 1: Percentage of ground states with angular momentum L in the IBM with random one- and 
two-body interactions obtained for N = 16 and 1,000 runs. 
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Figure 2: Probability distribution P(R) of the energy ratio R of Eq. (|l|) in the IBM with random one- 
and two-body interactions obtained for N = 16 and 1,000 runs. 
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Figure 3: Percentage of ground states with angular momentum L = (solid line), L = 1 (dashed line) 
and L = N (dotted line) in the vibron model with random one- and two-body interactions obtained for 
10 < N < 20 and 100,000 runs. 



10 











2 
R 



Figure 4: Probability distribution P(R) of the energy ratio R of Eq. (|7|) in the vibron model with random 
one- and two-body interactions obtained for N = 20 and 100,000 runs. 
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Figure 5: As Fig. [§ but for iV = 21. 
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Figure 6: Percentage of ground states with angular momentum L = (solid line), L — 1 (dashed line) 
and L = N (dotted line) in the vibron model with random one- and two-body interactions obtained in a 
mean-field analysis for 10 < N < 20. 



12 



